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We present the geometric structure of radiative transfer and the symmetries associated with its 
known conserved quantities. Our geometrization exploits recent work in the literature that enables 
the Hamiltonian formulation of radiative transfer to be obtained from a phase space representation 
of electromagnetic theory. Cosphere bundle reduction yields the traditional description over the 
space of "positions and directions", and geometrical optics arises as a special case when energy is 
disregarded. It is also shown that, in idealized environments, radiative transfer is a Lie-Poisson 
system with the group of canonical transformations as configuration space and symmetry group. 

PACS numbers: 42.15.— i, 46.05.+b, 47.10.— g 



Radiative transfer describes the transport of electro- 
magnetic energy in macroscopic environments when po- 
larization effects are neglected [1 . The theory originates 
in work by Bouguer [2j |3] and Lambert [4] in the 17 th 
century, where the measurement of light intensity was 
first studied systematically, and was later extended, for 
example to include transport effects [SHE]. Despite these 
efforts, however, to this day radiative transfer remains 
phenomenological, largely without connection to the fun- 
damental theories of light, and without a suitable mod- 
ern mathematical formulation — despite its importance in 
a multitude of fields, such as medical imaging, remote 
sensing, computer graphics, atmospheric science, climate 
modelling, and astrophysics. 

In this letter, we explain the physical foundations of 
radiative transfer in media with varying refractive in- 
dex and we study its Hamiltonian structure and sym- 
metries using a modern mathematical formulation, see 
Fig. [I] Exploiting recent work in the literature, we de- 
rive radiative transfer from electromagnetic theory us- 
ing an ansatz that resembles semi-classical analysis. The 
Wigner transform is employed to lift Maxwell's equa- 
tions from configuration space Q C R 3 to phase space 
T*Q. Restricting the dynamics on T*Q to a non-zero 
energy level, one obtains at the short wavelength limit 
a transport equation for polarized light, and further ne- 
glecting polarization leads to radiative transfer theory in 
a Hamiltonian formulation. This derivation uncovers the 
phase space light energy density i G Den(T*Q) as cen- 
tral quantity of radiative transfer theory, and shows that 
radiance, which classically plays this role, is meaningful 
only in the context of measurements, in the setting Lam- 
bert considered when he introduced the concept. With 
the Hamiltonian formulation of radiative transfer on 6- 
dimensional phase space T*Q, the classical 5-dimensional 
description is obtained by reduction to the cosphere bun- 
dle S*Q = (T*Q \{0})/R+ where time evolution is de- 
scribed by contact dynamics. The conserved quantity 
associated with this symmetry is the light frequency, un- 



covering the structure of a long known conservation law. 
Fermat's principle is obtained from the Hamiltonian for- 
mulation of radiative transfer through a non-canonical 
Legendre transform. Geometrical optics hence arises as 
a special case of the theory. When for a fixed environ- 
ment Q the light energy density t t G Den(T*Q) at time t 
is considered as a configuration of radiative transfer, and 
when one assumes that the Hamiltonian flow of the sys- 
tem is defined globally, the configuration space of radia- 
tive transfer becomes the group DifT can (T*(3) of canon- 
ical transformation. This makes ideal light transport a 
Lie-Poisson system, and reveals a remarkable similarity 
to ideal fluid dynamics, cf. Table [IJ The symmetry asso- 
ciated with the Lie-Poisson structure is the conservation 
of light energy density along trajectories in phase space, 
providing a modern rationale for the classical law of "con- 
servation of radiance along a ray" whose symmetry was 
previously not known. The geometrization and physical 
foundations of radiative transfer presented in this letter 
embed the theory into modern physics and they enable — 
and warrant — a fresh look at new and old questions. 

The physical foundation of radiative transfer theory is 
provided by the macroscopic Maxwell's equations, Hamil- 
ton's equations for electromagnetic field theory [9 , p. 24], 
which in a source free region Q C R 3 are given by 
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where e : Q —> R and /i : Q — > R are the electric per- 
mittivity and magnetic permeability, respectively. With 
U = (E,H) T ', Eq. Ill can be written as 

(js = ps jjs ( 2 ) 

where P £ is the Maxwell operator. In Eq. [2] we also in- 
troduced the scale parameter e = X/d n with A being the 
wavelength and d n the average distance over which the 
refractive index n = y/eji : Q — » R varies. Macroscopic 
environments are characterized by A <C d n , which war- 
rants considering the short wavelength limit e — » 0. The 
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FIG. 1. The physical foundations and geometric structure of radiative transfer theory. 



observable of interest for radiative transfer theory is the 
electromagnetic energy density 



£ £ (q,t) = \\U\\l^ 
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for e — >> and when polarization is neglected. The clas- 
sical approach to study short wavelength asymptotics 
is the Wenzel-Kramers-Brillouin (WKB) approximation. 
However, the ansatz yields well-defined solutions only un- 
til caustics form, when the approximation becomes multi- 
valued, and the initial conditions must satisfy the WKB 
form u £ (q,t) = a(q, t) e^fa,*)/^ As shown recently [TO] , 
these limitations can be circumvented by studying the 
short wavelength limit on phase space T*Q = R 3 x Q C 
R 3 . The electromagnetic field U £ can be lifted to T*Q 
using the Wigner transform, yielding a 6 x 6 "matrix" 
density W £ = W[U £ )(q,p). Eq. [2]then becomes 

W £ = -{{p £ ,W £ }} MB (4) 

where p £ = po + ep\ + e 2 p2 + ■ • • is the symbol of the 
Maxwell operator P £ , the dispersion matrix for the sys- 
tem [321 Eq. 3.63], and the bracket {{ , }} MB on the right 
hand side is the matrix-valued Moyal bracket [15j Chap- 
ter 16.3]. In contrast to the scalar case where the Moyal 

bracket satisfies {{ , }} MB £ ^°> { , } and reduces to the 
Poisson bracket { , } as e — » 0, in the matrix valued case 
one has 

{{p £ , W £ }} = l - (p £ W £ - W £ p £ ) + ■■■ + 0(e) (5) 

with the first term diverging as e —> since matrix mul- 
tiplication does not commute. For the short wavelength 
limit, this can be avoided by restricting the dynamics to 
an eigenspace, or allowed energy level, of the zero th order 
Maxwell symbol po? the analogue of the Bohr-Sommerfeld 
quantization condition [16 . The eigenvalues ofpo are [12] 
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with each having multiplicity two. Projection onto 
the eigenspace associated with a non-trivial eigenvalues 
A a , a G {1,2} and taking the limit e — >• 0, one obtains [11] 



W° a +{\ a ,W a }=[W° a ,F°] 



(7) 



where the projection W® of the limit Wigner distribution 
W° onto the 2-dimensional eigenspace is given by 



w° = n a w°n a 
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W® is formed by the Stokes parameters I,Q,U,V for po- 
larized light for the plane wave propagating in "direction" 
p G T*Q : and this also provides the physical interpreta- 
tion of p as the wave vector p = n(q)/cus, with uj being 
the angular frequency and s a unit vector. With Eq. [8j 
the right hand side of Eq.[7]can be interpreted as describ- 
ing the change in polarization during the transport by the 
Poisson bracket {A a ,W^}. Radiative transfer does not 
consider polarization. With the Stokes parameters for 
unpolarized light, Q,U,V = 0, Eq. [7] becomes 



e,H}. 



(9) 



The Hamiltonian H £ F(T*Q) in Eq. [9] is given by the 
eigenvalue A a , that is 



H(q,p)=± 



n(q) 



\\P\\ 



(10) 



so that the associated Hamiltonian flow is determined 
by the refractive index n(q). The trace £ = tr (I^ ) = 
£(q,p) dqdp G Den(T*Q) is related to the electromag- 
netic energy density S(q) by the fiber integral 



£(q)= f t=\ 
Jt*q Jt*c 



(ii) 



and the phase space density £ is hence the central quan- 
tity in radiative transfer theory, representing the en- 
ergy of the electromagnetic field at the short wavelength 
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limit when polarization effects are neglected. We denote 
£ G Den(T*Q) as phase space light energy density or light 
energy density. 

To relate £ G Den(T*Q) to radiance, the central quan- 
tity in the classical formulation of radiative transfer, we 
have to consider measurements, the question Lambert 
was studying when he introduced the concept [4] . Using 
standard results from tensor calculus, one obtains for the 
energy E flowing through a 2-dimensional surface hA in 
a time interval [ti,^] that [171 Chapter 3.2.6] 
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(12a) 



= fV ~T\ I C(q,p)(p-n)dAdpdt (12b) 

Jtt JM n W J T-M 

where p is a unit vector, n the surface normal of .M, 
and T~ M. the positive half-space of T*Q as defined by 
the surface normal n\ q . Parametrizing the light energy 
density in spherical coordinates [T7J Chapter 3.2.5], in- 
finitesimally the measurement can be written as 

C(q,p, v) (p ■ ft) dAdp = n-(C(q,p, v)pdAdpdv) . (13) 

When no measurement surface, and hence no normal n, 
is fixed an infinitesimal measurement is thus given by 



A = C(q,p, v) dA± dpdv 



where dA± (p) = p dA is the standard area form for a sur- 
face orthogonal to the flow direction p. The differential 
2-form A G fl 2 (Q) provides a modern interpretation of 
classical radiance, with the cosine term (p • n), prevalent 
in the classical literature, arising from the pullback of A 
onto a surface with normal n. 

Remark 1. It has often been overlooked that radiance 
is meaningful only in the context of measurements, which 
led to considerable confusion as in [18] , while the quantity 
that is naturally transported in radiative transfer is the 
phase space light energy density £ G Den(T*Q). 

Eq. [9] describes radiative transfer theory as a Hamilto- 
nian system on 6-dimensional phase space T*Q. In the 
literature, however, the theory is usually defined over the 
5-dimensional space of "positions and directions". The 
connection between the two descriptions can be obtained 
by considering the symmetry associated with the well 
known conservation of frequency. Since the momentum 
p is proportional to frequency the associated Lie group 
action is m a (q,p) = (q,ap) for the multiplicative group 
(R + ,-) acting on the fibers T*Q. Moreover, since the 
Hamiltonian in Eq. [10] is homogeneous of degree one in 
the momentum, the dynamics can be reduced from T*Q 
to the 5-dimensional cosphere bundle [I9j App. 4] 



S*Q = (T*Q\{0})/R + 



(15) 



where the time evolution is a contact Hamiltonian flow 
along the Reeb vector field [20] . The cosphere bundle 
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FIG. 2. Non-canonical Legendre transform relating radiative 
transfer to geometrical optics given by Fermat's principle. 

S*Q, which can be identified with the sphere bundle 
S 2 Q using a metric, provides a modern interpretation for 
the classical space of "positions and directions" , and the 
Hamiltonian being homogeneous of degree one explains 
why a description on S*Q is possible. 

Remark 2. A question open to this date is the relation- 
ship between radiative transfer and geometrical optics. 
The connection can be establied with the Legendre trans- 
form of Fermat's principle, cf. Fig. [2] As is well known, 
directly performing the transform for the Lagrangian 
L = n(q) leads to a vanishing Hamiltonian [2T] , Fol- 
lowing Arnold [T9] , we hence consider instead of length, 
given by L = n(g), the geometrical energy of a light path, 
given by L = n 2 (q). This La grangian can be interpreted 
(14) as a diagonal metric gij = n 2 (q)/c 2 5ij and the associated 



geodesic flow is equivalently described by the Hamilto- 
nian H = g l i pipj [22] p. 51] where p is the canonical 
momentum which is related to the kinetic momentum by 
p = ujp. Reverting the transition from path length to 
path energy, and including the factor of uo corresponding 
to energy we obtain for the Hamiltonian again Eq. [To] If 
we trace the flow in Fig. [2] backwards, we see that, from 
the point of view of radiative transfer, geometrical optics 
is a special case when energy transport is disregarded. 

A central result in classical radiative transfer theory 
is the "conservation of radiance along a ray" [23 . The 
symmetry associated with this conservation law becomes 
apparent when the light energy density £ t at time £, for 
all (q,p) G T*Q, is considered as one configuration of 
the system. Time evolution can then be described by 
the pullback £ t = n^io by the time t flow map rj t along 
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the Hamiltonian vector field Xh arising from Eq. 
When Xh is defined globally, the set of all maps n t forms 
the infinite dimensional Lie group Diff can (T*Q) [24] of 
canonical transformations, and with respect to £o all 
physically valid configurations £ t can be described by 
an element rj t in the group. Diff can (T*Q) can hence be 
considered as the configuration space of ideal radiative 
transfer, equipping the system with a Lie-Poisson struc- 
ture, cf. Fig. [3] The Lie algebra $ of Diff can (T*(5) can 
be identified with the generating Hamiltonian functions, 
= ^(T*^), in which case the dual Lie algebra 0* be- 
comes Den(T*Q) [26 . With $* ^ Den(T*Q), it is natu- 
ral to consider the light energy density £ as an element in 
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FIG. 3. The structure of Lie-Poisson systems. Classical exam- 
ples of such systems are the rigid body, where the Lie group is 
SO (3) (with a left invariant Hamiltonian), and the ideal Eu- 
ler fluid, where the Lie-group is the group Diff /X (Q) of volume 
preserving diffeomorphisms 24, 25] . 



0*. Time evolution is then described by coadjoint action 
Diff can (T*(5) x — » in the Eulerian representation, 
which infinitesimally is [T7J Sec. 3.3] 



t: 



31 



ad^i 



e,H}. 



(16) 



That the Poisson bracket describes infinitesimal coad- 
joint action ad* : q x — >> provides a rigorous justifi- 
cation for considering Diff can (T*Q) as the configuration 
space for ideal radiative transfer [27j Sec. 6]. In Eq, 
T~L = Ull] is the field Hamiltonian 
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H[£] 



I 

Jt*c 



t(q,P)H(q,p) dqdp 



(17) 



which is defined through the "single particle" Hamilto- 
nian H(q,p) in Eq. 
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With the light energy density 
as an element in the dual Lie algebra, it follows imme- 
diately from the general theory of Lie-Poisson systems 
that the momentum map Jr is the convective light en- 
ergy density and that it is conserved [9j Theorem 11.4.1], 
cf. Fig. [3] By the change of variables theorem, this 
can be interpreted as conservation of light energy density 
along trajectories in phase space, providing a modern for- 
mulation and justification for "conservation of radiance 
along a ray". Remarkably, with the Lie-Poisson struc- 
ture a close analogy between ideal radiative transfer and 
ideal fluid dynamics exists, cf. Table [l[ Next to the flow 
on phase space T*Q, one can obtain a functional ana- 
lytic description, related to the operator formulation of 
radiative transfer [28 , using Stone's theorem, yielding 
£ t = 77*^0 = Mt £0 where U t is a unitary operator. 

Our geometric formulation of radiative transfer and the 
identification of the Lie group structures underlying the 
known conservation laws clarifies and unifies earlier work 
in the literature. The use of tensor calculus overcomes 
the limitations of the current formulation, and also of 
earlier attempts that employed vector calculus [29] and 
measure theory [30] to obtain a modern mathematical 
foundation for radiative transfer. Our derivation of the 
Hamiltonian formulation of radiative transfer follows re- 
cent work in applied mathematics [TTJ [12] that employs 
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TABLE I. Correspondence between ideal fluid dynamics and 
ideal radiative transfer. The fluid velocity is denoted by v G 
3£div(Q) and uj G ^ 2 (Q) is the fluid vorticity. 

the Wigner transform to study short wavelength asymp- 
totics, chiefly because of the simple and linear relation- 
ship in Eq. [TT] The use of the Wigner transform parallels 
earlier approaches [31] [32] . The authors of [TTJ [12] put 
these on a rigorous mathematical basis, which is crucial 
to treat the delicate questions arising from the highly os- 
cillatory nature of Z7 £ , and they elucidated many aspects 
that previously seemed serendipitous. Currently lacking 
is a treatment of the transition from Maxwell's equa- 
tions to radiative transfer from an algebraic and geomet- 
ric point of view. Gerard and coworkers [TT] . for example, 
did not interpret their equations as matrix- valued Moyal 
bracket, and hence did not relate their results to the large 
body of work on quantization [35] . The status and short- 
comings of many classical derivations of radiative transfer 
theory were recently summarized by Mishchenko [38] . 

An interesting extension of our geometric formulation 
would be the inclusion of dissipation and scattering. The 
classical Beer-Lambert law is naturally compatible with 
our description [TTJ Remark 3.29] but it would be inter- 
esting, and of considerable practical relevance, to extend 
this to more general scattering models, for example based 
on recent work on kinetic theory [39]. An important ob- 
jective for future work should also be to develop the con- 
nection between our formulation of radiative transfer and 
work by Mishchenko [40 , 41 where radiative transfer the- 
ory for media with scattering particles was derived from 
electromagnetic theory using an alternative approach. 
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